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ABSTRACT 

We present results for the assembly and star formation histories of massive L*) red sequence 
galaxies in 11 spectroscopically confirmed, infrared-selected galaxy clusters at 1.0 < z < 1.5, the 
precursors to present-day massive clusters with M ~ lO^^M©. Using rest- frame optical photometry, 
we investigate evolution in the color and scatter of the red sequence galaxy population, comparing 
with models of possible star formation histories. In contrast to studies of central cluster galaxies at 
lower redshift (z < 1), these data are clearly inconsistent with the continued evolution of stars formed 
and assembled primarily at a single, much-earlier time. Specifically, we find that the colors of massive 
cluster galaxies aX z k, 1.5 imply that the bulk of star formation occurred at z ~ 3, whereas by 
z K, 1 their colors imply formation at z ~ 2; therefore these galaxies exhibit approximately the same 
luminosity- weighted stellar age at 1 < z < 1.5. This likely reflects star formation that occurs over an 
extended period, the effects of significant progenitor bias, or both. Our results generally indicate that 
massive cluster galaxy populations began forming a significant mass of stars at z > 4, contained some 
red spheroids by z « 1.5, and were actively assembling much of their final mass during 1 < z < 2 
in the form of younger stars. Qualitatively, the slopes of the cluster color-magnitude relations are 
consistent with no significant evolution relative to local clusters. 

Subject headings: galaxies: clusters: general - galaxies: elliptical and Lenticular, cD - galaxies: evolu- 
tion ~ galaxies: formation - galaxies: photometry 



1. INTRODUCTION 

A key goal in observational cosmology is to measure 
the formation and assembly of cosmic structures, a his- 
tory that can be traced by examining the stars in massive 
galaxies. The simplest model consistent with many ob- 
servations is that the stars in massive early-type galaxies 
(ETGs) formed in a single short burst at high redshift 
and evolved passively (i.e., with no further star forma- 
tion) thereafter. The dramatic expansion of data from 
high-quality wide-area galaxy surveys and deep lookback 
studies has permitted the rigorous testing of this pic- 
ture, revealing a more complex formation and assembly 
history. 

Present-day clusters contain a tight red sequence of 
galaxies, indicating early a nd uniform star formation 
histories (jBower et al.lllQQ^ . As fossil records of clus- 
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ter assembly, red sequences were found t o evolve nearly 
passively since z ^ 1 ( Aragon-Salaman ca et al.l 119931 : 
[Stanford et aLlfiggllM^zin et al . 2008), suggesting neg- 
ligible ongoing star formation. Their colors imply for- 
mation at z ^ 2 to 3, and studies at z ^ 1 (e.g., 
Blakeslee etld] 120031: iMei et al.l l2006bl: iLidman et alj 



2008HMei elall [20091 : IStrazzullo et al.ll2010t ) found con- 
sistency with this picture. The single short burst 
model predicts that the red sequences become bluer and 
wider nearer their form ation epoch, a t r end f ound by 
iHilton et al] (|2009f ) and iPapovich et"all (|2010i) in two 
clusters at z « 1.5 — 1.6. 

This passive view is brought into que stion by some ob- 
servat ions of clusters at z > 1. While [Eisenhardt et al.l 
([2008) found that the colors of cluster galaxies at z < 1 
are consistent with the passive evolution of a group of 
stars formed in a short burst at Z f = 3, the z > 1 candi- 
dates favor an earlier formation. iMancone et al.l (|20100 
found that cluster galaxies' characteristic 3.6/xm and 
4.5/im magnitudes (from the Spitzer Deep, Wide-Field 
Survey (SDWFS, .Ashbv et al . 2009,)) evolve in a manner 
consistent with little mass assembly at z < 1.3 (see also , 
[Lin et al.|[2006llMuzzin"eraI|[2008l[Rettura et al.|[20Tl[) . 
but are systematically fainter than th e extrapolation of 
passi ve models at z > 1.3 (see also, iFassbender et al.l 

IMB. 

Furthermore, wc now know the red sequence is a 
snapshot of the currently quiescent galaxies, and that 
galaxies evolve both o n and off during formative events. 
Such transform ations (" Butcher fc Oemieilll978l : iDressleil 
[19801 iTran et a l. 2003; Mcintosh et al.ll2008D may cause 
the red sequence's average stellar age to evolve more 
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Table 1 

Summary of HST observations of JSCS clusters. 



INamc 


T5 A 
IX. A. 


JJec. 


z 


L/iViJJ Area 


Instrument / r liter 


JJeptn 


JJeptn 










(arcmm ) 


(optical) 


(optical) 


( Wl Co/l" 160W) 


ISCS J1429. 2+3357 


14:29:15.16 


33:57:08.5 


1.059 


3.1 


WFPC2/F814W 


26.1 


24.7 


ISCS J1432.4+3332b 


14:32:29.18 


33:32:36.0 


1.112 


3.6 


ACS/F775W 


26.2 


24.8 


ISCS J1426.1+3403 


14:26:09.51 


34:03:41.1 


1.136 


3.3 


WFPC2/F814W 


26.1 


24.7 


ISCS J1429. 2+3425''''^ 


14:29:15.16 


34:25:46.4 


1.162 


1.7 


WFPC2/F814W 


26.6 


24.7 


ISCS J1426.5+3339 


14:26:30.42 


33:39:33.2 


1.163 


1.9 


WFPC2/F814W 


26.4 


24.5 


ISCS J1427.9+3430'= 


14:27:54.88 


34:30:16.3 


1.235 


3.3 


WFPC2/F814W 


26.7 


24.7 


ISCS J1434.5+3427 


14:34:30.44 


34:27:12.3 


1.243 


4.4 


ACS/F775W 


26.1 


24.5 


ISCS J1429.3+3437 


14:29:18.51 


34:37:25.8 


1.262 


4.2 


ACS/F850LP 


26.3 


24.4 


ISCS J1432.6+3436 


14:32:38.38 


34:36:49.0 


1.349 


4.4 


ACS/F850LP 


26.0 


24.8 


ISCS J1433.8+3325 


14:33:51.13 


33:25:51.1 


1.369 


4.5 


ACS/F850LP 


26.1 


24.8 


ISCS J1434.7+3519 


14:34:46.33 


35:19:33.5 


1.372 


4.5 


ACS/F850LP 


26.1 


24.8 


ISCS J1438.1+3414 


14:38:08.71 


34:14:19.2 


1.413 


4.5 


ACS/F850LP 


26.1 


24.7 


ISCS J1432.4+3250 


14:32:24.16 


32:50:03.7 


1.487 


8.0 


ACS/F814W 


26.0 


25.8 



^ 5-sigma point source depth in 0.8" diameter; AB system 

^ HST pointing did not intersect cluster core due to guide star requirements 



No single identifiable CMR 

slowly (jvan Dokkum fc Franxl I2001D than predicted by 
the single-burst model. Therefore, estimates of stellar 
age from the red sequence color and color dispersion 
(|Bower et al.lll99^ iDemarco et al.ll2010l ) can provide im- 
portant constraints on galaxy assembly in clusters. 

In this work, we report Hubble Space Telescope 
[HST) follow-up imaging of cluste rs drawn from the 
infrared-selected sample of Eiscnhar dt et al.l (j2008D . the 
Spitzer/mkC Shallow Cluster Survey (ISCS). This sur- 
vey used near-IR imaging, an excellent proxy for stel- 
lar mass, to select structur es that are the prec ursors 
to today's massive clusters (jBrodwin et al.|[2007D . We 
directly study their histories by examining the red se- 
quence galaxies in 1 1 spectroscopically confirmed clusters 
at 1.0 < z < 1.5. 

In ij2] we describe the candidate selection and data 
collection, and in fj3] we quantify the rest-frame optical 
color-magnitude relations. We construct simple models 
for their evolution in fj^l and we compare these models 
with our measements in fj5] We discuss the implications 
of these results in iJB]and conclude in iJTl Throughout this 
work our observed-frame magnitudes are on the AB sys- 
tem, but for convenience in comparing to previous work 
we also calculate rest-frame U-V colors on the Vega sys- 
tem. We assume a flat cosmology with Oa = 0.72 and 
h = 0.71, consistent with ILarson et al.l (|2011l . WMAP7). 

2. DATA, REDUCTION, AND GALAXY CATALOGS 

2.1. Cluster Catalog 

The Syitzer/mkC S hallow Cluster Survey (ISCS; 
lEisenhardt et al.l |2008() identified cluster candidates 
spanning 0.1 < z < 2 over a 7.25 deg^ area in 
the Bootes field of the N OAO Deep Wide-Field Sur- 
vey pannuzi fc Devi |1999() . The clusters were identi- 
fied as three-dimensional overde nsities using accurate 
optical/IR photometric redshifts (jBrodwin et al.l [2006[ ) 
calculated for the 4.5 /xm flux-l imited (8.8 /LtJ y at 5cr) 
IRAC Shallow Survey (Eiscnh ardt et al.l [20M) catalog. 
Thus, the clusters are selected primarily by stellar mass, 
largely independent of the presence of a red sequence. 
There are 335 ISCS cluster and group candidates, 106 
of which are at z > 1. Over 20 clusters at z > 1 have 
now been spectroscopically confirmed with W. M. Keck 



Telescope and/or HST spectroscopv (IStanford et al"] 
'2005', 2012: Elston ct al.' 20CI6t fBrodwin et al.ll2006ll201ir 
lEisenhardt et al..,2008. Zeimann et al. in prep., Brodwin 
et al. in prep.). 

2.2. Imaging and Reduction 

In this work we present sensitive, high-resolution HST 
follow-up imaging obtained from a variety of programs 
(GO proposal IDs 10496, 11002, 11597, and 11663) for 
a subset of the confirmed z > 1 ISCS clusters. In or- 
der to construct a rest-frame optical color magnitude 
diagram with the most leverage on the galaxies' stel- 
lar ages, instrument /filter combinations were chosen to 
bracket the 4000A break. Specifically, each cluster was 
observed in th e near-infrare d with the Wide Field Cam- 
era 3 fWFC3: lKimble et al.lf2008) filter F160W. Optical 
data w ere taken with th e Advanced Camera for Surveys 
(ACS; IFord et al.l [19981) in filters F775W and F850LP 
or F814W, or with the Wide Fie ld Planetary Camera 
2 (WFPC2: iHoltzman et aL|[T995[) filter F814W, ah ap- 
proximating broad / band or I+z band coverage. We 
provide a summary in Table [T] 

The HST images were reduced using stan dard proce- 
dures with the MultiDrizz le software (Koek emoer et afl 
[2OOI iFruchter et all 120091) . The WFC3/F160W imag- 
ing was obtained in a single pointing of 4xl03s (dithered 
between each exposure) for all but one of the clusters. 
For cluster ISCS J1432.4-h3250 (z = 1.487) there were 
two slightly overlapping pointings, each consisting of 700s 
of integration time in a two-dithered pattern. The na- 
tive (undersampled) WFC3 0.13" pixel scale was driz- 
zled down to 0.065" to better match the pixel scale of 
ACS/WFPC2 and to better recover the point-spread 
function. 

The regions from which we construct CMDs are limited 
to the overlap between the optical and WFC3/F160W 
fields, 2 to 8 arcmin^, corresponding to fields of view 
^ 0.7 to 1.5 Mpc in diameter. Therefore, all results that 
follow pertain only to the cluster cores, and we will re- 
fer to galaxies within ~ 1 Mpc as "central". The HST 
fields for these clusters differ in total size, and so we 
have checked that the conclusions of this paper are un- 
changed when we consider only sources within a fixed 
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0.7 projected physical Mpc radius from the spectroscopic 
cluster centers. 

The ACS/F814W observation of ISCS J1432.4+3250 
(z = 1.487) is comprised of 8x564s exposures. The 
ACS/F850LP and F775W images for seven clusters 
were obtained as part of the HS T Cluster Supernova 
Survey (PI Perlmutter, GO-10496 rPawson et all l200l 
iSuzuki et all [Ml iBarbarv et all l2012allbf ). The pro- 
gram consisted of 219 orbits targeting rich clusters at 
z > 1, and the images comprise multiple overlap- 
ping pointings covering tens of arcmin^. The single- 
pointing WFPC2/F814W images were obtained in a 
joint Spitzer / HST program to survey Luminous Infrared 
Galaxies (LIRGs) in clusters at 2; > 1 {HST proposal ID 
11002, Spitzer proposal ID 30950). 

Scamp (Bcrtin 2006) was applied to align the world co- 
ordinate system (WCS) of each image with the reference 
frame of SDWFS (itself matched to 2MASS). SWarp 
()Bertin et al.ll2002[ ) was used to resample the HST opti- 
cal images to match the pixel scale of the drizzled WFC3 
images and to project the images to the tangent plane. 
This process aligned the WFC3 and optical images so 
that pixels in all images for a given cluster covered pre- 
cisely the same regions of the sky. 

2.3. Photometric Catalogs 

From the reduced and registered images, we mea- 
sure fixed-aperture photome try using Source Extractor 
(SE; iBertin fc Arnout"sl 11996') in dual-image mode with 
sources detected in the WFC3/F160W images. For com- 
parison we compute colors using aperture diameters of 
0.4" to 1.5", and choose 0.8" for our default color as 
a reasonable balance between preventing interference in 
our crowded fields and intersecting a physically represen- 
tative portion of each galaxy. We also constructed colors 
using galaxy-specific apertures scaled to the half-light ra- 
dius. We manually inspected the CMRs made with each 
aperture, and found little change in the appearance and 
location of the red sequences. Our results are thus in- 
senstive to our choice of photometric aperture. We use 
SE's MAG.AUTO in the F160W band as an approximate 
total magnitude for our CMRs. 

Differences in ETG color profiles may cause a bias in 
CMR properties owing to our fixed-aperture color pho- 
tometry. The aperture tests described above suggest that 
the intrinsic widths of the color distributions dominate 
this bias for most galaxies. Furthermore, in what fol- 
lows we assume that the slope of the rest-frame CMR 
matches th at of the Coma cluster, in which the slope was 
measured (jEisenhardt et al.ll2007D in an aperture corre- 
sponding to a nearly identical physical scale as the one 
used herein. If high-redshift ellipticals form the cores of 
present-d ay ETGs (e.g . , as suggested by the measure- 
ments of iTruiillo et al.l 120111 ). then we expect minimal 
bias owing to this effect. 

Our photometric uncertainties are dominated by sky 
shot noise, which we estimate in each image from the dis- 
tribution of fluxes through 5000 randomly-placed aper- 
tures with a diameter matched to our catalogs. The flux 
distribution is to a good approximation normal except for 
the sources contaminating the positive half. We compute 
the background level dsky by fitting a normal distribu- 
tion to the negative half. We verified this procedure by 
confirming that it produces the correct scaling in pho- 



tometric scatter of sources detected in the sets of dither 
images before making the final stack. 

We calculated the standard completeness of the 
WFC3/F160W-selected catalogs by generating artificial 
objects with known total magnitude and i?^/^ light pro- 
files where Re is distributed as a Rayleigh function with 
aR^ = 3.5 pixels (~ 0.23 arcseconds, or ~ 2 Kpc at 
z — 1.5), matching well the observed SE-derived Re 
distribution for sources in these fields. We found that 
the WFC3/F160W-detected catalogs for each cluster are 
more than 90% complete at an AB magnitude of 23.5, 
and every field is roughly 97% complete at H*{z) + 1.5, 
which is approximately H = 22-23 for these redshifts. 
H*{z) -1-1.5 is the magnitude cut we impose for the CMR 
analysis presented in this work (Section 13. 2|) . 

3. COLOR-MAGNITUDE RELATIONS 

We aim to study the star formation epoch of ETGs on 
the red sequences for our z > 1 cluster sample using the 
color and color scatter of the red sequence CMRs, which 
are identified using morphological classification and color 
selection techniques. The color and scatter of a CMR 
serve as proxies of the formation history of the galaxies 
that comprise it. The simplest physical parameters that 
they probe is the average formation time or age of the 
galaxies' stars. If the galaxies formed separately in short 
bursts over a period of time, then the color scatter will 
be inversely proportional to the time since they stopped 
forming, whereas if they formed in a synchronized fash- 
ion, then the scatter between galaxies may be always 
small. Since this is true regardless of their average age 
(and hence average color) , then we are motivated to mea- 
sure both color and scatter in order to better constrain 
the history of these galaxies. 

3.1. Morphological Classification 

Visual inspection of the individual galaxies in the 
F160W-selected catalogs was carried out by one of us 
(S. A. S.) using the WF C3 images. T-types (RC3; 
Ide Vaucouleurs et al.lll99l[ ) were assigned where possible 
to all objects down to FIQQW = 23.5. If the object was 
too faint to accurately classify, too close to a bright star, 
or affected by bad data, no morphology was assigned. In 
all other cases an attempt was made to determine the T- 
type on the scale ranging from —5 for a giant elliptical, 
to 10 for an irregular galaxy. In addition, information 
was recorded on the environment of the object, if it was 
compact or appeared inclined to the line of sight, and if it 
appeared to be interacting with any neighboring objects. 

Quantitative measurements of galaxy shapes were 
made by fitting a single Sersic profile to every ob- 
ject in the F160W- selected catalogs . The fitting was 
done with Galfit jPeng et al.ll2010( ). and Galapagos 
(jHaufiler et al.l[20Tll ) was employed to automate the task 
of running Galfit. Galfit does two dimensional fitting 
and can fit an arbitrary number of models to any image 
while accounting for the image's point spread function. 
Galapagos automates the task of running Galfit by 
generating first-guess Sersic parameters from an SE run, 
measuring the sky for each galaxy, deciding when to si- 
multaneously fit nearby galaxies, and deciding when to 
subtract off brighter objects that might bias a fit. It then 
uses Galfit to fit a single Sersic profile to every object 
in the image and builds a catalog of the results. 
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Figure 1. Comparison between morpiiological selection techniques for clusters ISCS J1429. 2+3357 (z = 1.059) and ISCS J1438. 1 +341 4 
(z = 1.413). We identify objects in the top and bottom panel based on visually-classified and a utomated morphologies, respectively ( ^3.111 . 
The black dotted lines show the CMR calculated from the Coma cluster I IEisenhardt et al]|2007l ) and evolved (assuming a single short burst 
at Zf = 3) t o the observed redshift. We plot a vertical dotted line a,t H = H* + 1.5, representing the magnitude selection limit we apply in 
Section |3 21 Gray crosses represent all spectroscopically confirmed cluster members currently known that are covered by the HST images 
and that fall within the color-magnitude ranges of the CMD. An updated description of the spectroscopy will be presented by Zeimann 
et al. (in prep.) and Brodwin et al. (in prep.). Where visible in the upper right corners, the diagonal blue dashed lines approximate 
the boundary of the 5(t color and F160W magnitude measurements. For ISCS J1438. 1+3414 {z = 1.413), we note the presence of many 
early-type galaxies ^ 1 magnitude blueward of the identified cluster red sequence. Several of these objects are part of a foreground structure 
evident in the photometric redshifts of the ISCS catalog, and are offset from the spectroscopic cluster center by ~ 1 arcminute. Therefore 
we believe those relatively blue early-type galaxies are interlopers. We find that both morphology classification schemes yield similar CMRs 
for all of our clusters, so we present the CMDs in FigurefJlEl and[4]using only the visual typing, and analyze them in Figure|4] 
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Figure 2. CMDs for our cluster sample, in ascending redshift order. These graphs follow th e sam e formatting as the top panels of the 
two CMDs in Figure [T] where plot symbols represent the visual classification system of Section [XT] 



In Figure [T] we present a comparison of the visual 
and quantitative morphologies in CMDs of two clusters 
that span the redshift range of our sample. The early- 
type CMRs resulting from these methods are qualita- 
tively quite similar in these two clusters, as well as in 
the rest of our sample, and lead to quantitatively similar 
results. For clarity we therefore present the remainder 
of our CMDs using only the visual T-types, but for com- 
pleteness retain a comparison with the quantitive mor- 
phologies in Figures [1] and [6] and in several tables. 

3.2. Red Sequence Selection 

We use the colors, magnitudes, morphological classifi- 
cations, and available spectroscopy to identify early-type 
members of the red sequence in each cluster. Multi-band 



photometric redshifts, which are not considered in this 
work, are available for only a small fraction of the galax- 
ies selected in the WFC3/F160W images, because these 
-ffS'T-based catalogs are dominated by fainter galaxies in 
fields too crowded for reliable ground-based photometry. 

The basic observables, the CMDs for 13 clusters, are 
presented in Figures [2] and [3l In each panel we plot 
the Coma CMR ( Eisenhard t et all 120071 ) evolved to the 
appropriate redshift (see Section |4] for details) . 

In ISCS J1429.2+3425 (z = 1.162) and 
ISCS J1427.9-H3430 {z = 1.235), shown in Fig- 
ure [3l there exist multiple known structures at different 
redshifts in the HST field of view. Moreover, the 
intersection between the WFC3/F160W image of 
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Figure 2. Continued from previous page. 

ISCS J1429.2+3425 (z = 1.162) and its WFPC2/F814W 
image is small and does not include an identifiable 
CMR. To avoid a possible bias from including projected 
large scale structure, as opposed to massive, collapsed 
clusters, we omit these two clusters from our analysis. 
It is possible that one or two other members of the 11 
remaining clusters are similar outliers, but lacking such 
evidence at present, we cannot confidently omit them 
from our analysis. 

For the clusters with a clear red sequence, the rest- 
frame CMR slope is within 10% of that in Coma, while 
the color zeropoint and scatter clearly vary. The shal- 
low slope of the CMR has been found to be nearly con- 
stant with redshift (Kodania & Arimoto 1997), and thus 
has been interpreted as evidence for a uniformly old se- 



quence in mass and metallicity (e.g., |Larson||19740, and 
not one based on galaxy age. If the tilted CMR is entirely 
the result of a mass-metallicity relation with a similar 
slope froin z = 2 to (as imp lied by the results of e.g., 
lErb et all [20061: iBrooks et ani2007[) . then we do not ex- 
pect cosmic evolution in the CMR slope so long as we can 
satisfactorily select galaxies that satisfy this relation. In 
what follows we fix the CMR slope in each cluster's rest 
frame, treating it as independent of age, and consider 
only measurements of the color zeropoint and scatter. 
We will discuss possible deviations from this assumption 
as needed. 

We isolate CMRs in the 11 clusters using the following 
procedure. First, we subtract from every galaxy's color 
our fiducial Coma CMR model (a single stellar popula- 
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Figure 2. Continued from previous page. 

tion formed in an instantaneous burst at zj = 3). We 
plot the resulting quantity A, in magnitudes, in Figure|31 
For all clusters except ISCS J1426. 5+3339 (z = 1.163), 
we define a red sequence galaxy (RSG) as having and 
A^^ - 0.25 < A < 0.75, where A^^ is the color differ- 
ence between our Zf — 3 assumption and a single stellar 
population formed at redshift Zf. We use Zf = 2 to 
set this lower limit on A instead of z/ = 3 in order to 
successfully select red sequences that may have formed 
more recently. This A range is chosen as a reasonable 
balance between including member RSGs and exclud- 
ing blue interlopers or late-type members that are not 
on the cluster's red sequence. This range of A, within 
which we select the RSGs, is listed for each cluster in 
Table [5J in all cases this is > 1.2 magnitudes. Except for 



ISCS J1426.5-I-3339 (z = 1.163), which has a very blue 
red sequence, this range cleanly selects the RSGs, and 
in particular the one composed of ETGs, that are visu- 
ally identifiable in Figures [HI For ISCS J1426.5-K3339 
(z ~ 1.163), we have manually modified the color selec- 
tion region to identify the RS: A2 - 0.5 < A < 0.75. 
It is unclear why ISCS J1426.5-h3339 (z = 1.163) ap- 
pears so much younger than other clusters at similar 
redshifts, though we have verified that our conclusions 
are unchanged if we omit it from our analysis. 

Next, we define the cluster RSGs as galaxies brighter 
than H*{z) -f 1.5 and selected by the above color cri- 
teria, where H*{z) is calculated by evolving the char- 
acteristic brightness for Coma (de Propris et al. 199^ 
to the cluster's observed redshift, assuming the galax- 
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Figure 2. Continued from previous page. 

ies formed in a single short burst at Zf = 3. Some 
of the RSGs selected in this way are spectroscopically 
confirmed cluster members, and we plot these as gray 
stars in all CMDs. However, most do not have accom- 
panying spectroscopic or photometric redshift informa- 
tion, and so the RSG samples may suffer from inter- 
loper contamination. In addition to the sample of all 
RSGs, we analyze two samples drawn from them to iden- 
tify likely ETG s: visually-selected E's or SO's as defined 
in Section 13.11 and quantitatively-selected ETGs with 
ns > 2.5. The latter choice is reaso nably consistent wit h 
selecting bulge-dominated objects (iSimard et al.ll201lD. 
as well as visually-typed E and SO galaxies (jSimard et al.l 
12009). 

Finally, to limit the influence of interlopers that are ex- 
treme outliers from the cluster's red sequence, we drop 
galaxies from the final ETG or RSG catalog that are fur- 
ther than two median absolute deviations in color from 
their A zeropoint, as measured below. Table [2] lists the 
number of galaxies satisfying our red sequence selections, 
as well as the median color uncertainty (for the sample 
of all RSGs) in the observed bands. 

3.3. Color Zeropoints and Scatters 

For each of the three ETG or RSG samples de- 
fined above, we measure the color zeropoint and in- 
trinsic scatter, and their associated uncertainties, ap- 
plying the biweight est i mates of location and scale 
(|Mosteller fc TukevlfTOTTt iBeers et al.l [T990l ) for the ze- 
ropoint location Aq and scatter about this location (ta. 
These two measurements define our GMR for each sub- 



sample. We calculate the color Co at H*{z) by adding 
back to Ao the previously-subtracted fiducial Zf — ?> 
CMR model, and the intrinsic scatter aint by subtract- 
ing in quadrature the median photometric error from the 
biweight scale estimate ua of the sample in question. 

Uncertainties are derived from 1000 bootstrap resam- 
plings, from which we measure the scale of the resulting 
Co and (Tint distributions using the median absolute devi- 
ation (MAD). Specifically, ac = 5mad = MAD/0.6745, 
since a normally-distributed quantity will have MAD = 
0.6745a. Where the samples overlap, our measure- 
me nts of uu-y are c onsistent with those measured 
bv IMevers et all (l20T2h . except for ISCS J1438.1+3414 
(z = 1.413). At z > 1.3, the addition of the HST 
WFC3/F160W data permits an improved k-correction, 
and so we believe the present calculation to be reliable. 

For each cluster we present Co, uc, and their uncer- 
tainties in Table [31 The observed-band (Tint values are 
all significantly smaller than w 0.6 magnitudes, which is 
the value we would expect to measure if there were no 
overdensity (e.g., no red sequence) inside the A selection 
region. This implies that our measurements are driven 
by the real red sequences and are not an artifact of our 
color selection choices (Table [2]). 

4. SIMPLE MODELS OF COLOR AND SCATTER 

To associate the measured CMRs with formation 
epochs, we compute grids of m odel colors in our observed 
passbands using the models of iBruzual fc Chariot! (|2003l 
2007 version). Each model galaxy forms its stars in a 
single exponentially-decaying burst with a timescale of 
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Figure 3. CMDs of two clusters for which there are known structures at different but similar (2 ~ 1) redshifts. The symbols are the same 
as in Figure [2] In order to avoid potential confusion owing to multiple red sequences, these two clusters are not used in the subsequent 
analyses. 



Table 2 

Number of cluster red sequence galaxies. 









Number on 


OMR 


Spectroscopic Members'" 










Vis. 


Quant. 


All 


WFC3 


All 


median a„^c,t 


Name 


2 


A width'' 


ETGs 


ETGs 


RSGs 


Catalog 


RSGs 


All RSGs 


ISCS J1429. 2+3357 


1.059 


1.19 


13 


16 


26 


10 


5 


0.04 


ISCS J1432. 4+3332 


1.112 


1.25 


5 


7 


11 


6 


3 


0.05 


ISCS J1426. 1+3403 


1.136 


1.24 


15 


20 


30 


8 


6 


0.03 


ISCS J1426. 5+3339 


1.163 


1.76 


18 


12 


39 


6 


6 


0.04 


ISCS J1434.5+3427 


1.243 


1.25 


19 


19 


27 


10 


3 


0.05 


ISCS J1429. 3+3437 


1.262 


1.22 


23 


21 


33 


13 


8 


0.04 


ISCS J1432. 6+3436 


1.349 


1.24 


12 


16 


27 
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5 


0.05 


ISCS J1433. 8+3325 


1.369 


1.24 


9 


9 


24 


6 





0.07 


ISCS J1434.7+3519 
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1.24 


10 


14 


27 
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0.04 


ISCS J1438. 1+3414 


1.413 
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20 
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ISCS J1432. 4+3250 


1.487 


1.24 


15 
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0.10=^ 



Color width in magnitudes used to define the RSG selection region; see Section 13.21 

Spectroscopy for parts of this sample were described by [Stanford erall l(20M l: lEIston et all I^OOB ): IBrodwin eFan 1(2003. [200 ): 
lEis enhardt et aLl 1 I2OO81) . Updated membership identification, including those confirmed with the HST WFC3 grism, will be described 
by Zcimann et al. (in prep.) and Brodwin et al. (in prop.). 

Although this cluster has the deepest WFC3/F160W image in the sample, its red sequence galaxies are fainter (and redder), leading to 
a high average color uncertainty. 
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Figure 4. CMD residuals from Zf = 3 model and subsequent analysis for ea ch c luster, in ascending redshift order. We restrict attention 
to elliptical and SO galaxies satisfying the color and magnitude cuts of Section l3.2l All RSGs and the visually-selected ETG subsample are 
plotted as gray and black fille d circles, respectively. We plot a vertical dotted line at H = H* + 1.5, representing the magnitude selection 
limit we apply in Section 13.21 The location and scale measurements are shown as pairs of horizontal lines centered on the final measured 
Ao zeropoint value and spaced to surround the one-sigma region defined by the final measured (Tint scatter value. The blue solid lines 
correspond to the scatter of the subsample containing all RSGs, and the red solid lines correspond to scatter in the visual ETG subsample. 
Histograms a re ta bulated for each subsample on the right of each CMD, and include dashed black lines to show the color selection region 
set in Section 13.21 
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Figure 4. Continued from previous page. 

r = 0.1 Gyr and evolves passively thereafter until it is 
observed. We will label these models our "short-burst 
stellar population" (sbSP) models, whose evolution after 
a few hundred Myr is very similar to that of a simple 
stehar population (SSP). To construct the gri ds we use 
the EZGal code of lMancone fc Gonzaled (|2012l Fl. focus- 
ing on so l ar-me tallicity stellar populations drawn from a 
IChabrieil (|2003[ ) initial mass function (IMF), using for- 
mation epochs evenly spaced in time between Zf — 1.8 
and Zf = 10. 

For ease of comparison to other work, we convert our 
measu red quantit i es to th e rest-frame UBV filter set fol- 
lowing [MeTeE^D (|2006bD using these models. First, we 

http://www.baryons.org/czgal/ 



evolve a set of sbSP models with varying metallicity and 
formation redshift down to each cluster's redshift. Then 
we fit the relation {i — H)z — a + b{U — V)o to estimate 
the parameters a and b. The slope b is used to derive the 
rest-frame color scatter given an observed color scatter; 
specifically, auv — <yiH/b. We also use these relations to 
comp ute the slope of the Goma CMR (jEisenhardt et all 
120071) as it would be observed in the HST bands at these 
redshifts; these were used for selecting red sequences in 
Section 13.21 We visually inspected the linear fits and 
found good agreement in all cases. We also performed 
the analysis of Section [5] using the colors as-observed and 
using only directly comparable models at the appropriate 
redshifts. We found that our results are identical in both 
cases, confirming the validity of the conclusions based on 
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the k-corrected rest-frame colors. 

We construct model red sequ ence ETG color scat- 
ters followin g prev ious studies (iPapovich et al.l 120101 : 
iHilton et all l20Q9t iMei et al.l 120091 and references 
therein). We consider the simplest models in which the 
dispersion in optical ETG color for a single model set 
depends on the redshift at which star formation ceased, 
Zend- Prior to Zond, our models assume that cluster galax- 
ies form uniformly in time starting at zq = 10. In prac- 
tice, for each stellar IMF and metallicity we draw 10^ 
individual galaxy models with random formation times 



between zq and Zond and evolve them to each cluster's 
observed redshift. We then compute three quantities: 
the color dispersion and median color using t he biw eight 
scale and location estimators (iBeers et al.l Il990( ) , and 
the luminosity-weighted average formation time, which 
we label (z/), a typical (star-)formation redshift. The 
epochs Zend and (zy) describe these galaxies' last ma- 
jor star-forming episode, the nature and timing of which 
we seek to trace. In what follows we denote this set of 
models as our "constant SFR" (cSFR) models. 
To compare the color zeropoints between the observed 
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clusters and model galaxies, we measure the color of 
the observed red sequence at the characteristic bright- 
ness, H*(zohs)^ that w e evolve from the value in Coma 
(|de Propris et al.lll998D . We assume each model galaxy 
has the same mass, and thus refer to all CMR quanti- 
ties as those measured at H*{zohs)- A galaxy luminosity 
function ev olving by 0.5 magnitudes in the fashion of 
iMancone^ t al. (2010) will only mildly affect our inter- 
pretation of their stellar populations. Specifically, the 
comparison between observed and model color depends 
only weakly on the magnitude at which we compute nu- 
merical values of the colors from the derived CMRs. This 
owes to the CMR slopes in our observed bands being shal- 
low (w —0.1), so a 0.5 magnitude shift in the value of H* 
alters our derived U — V colors by only w 0.03 magni- 
tudes. To recover the results we obtain in Section [S] with 
a fixed formation epoch, H* must increase by > 3 mag- 
nitudes between z w 1.5 and z w 1, a change we believe 
is too large. We conclude that this effect may slightly 
bias our color measurements to the red but it does not 
change our qualitative results. 

The true star formation and assembly histories of 
ETGs, as observed in cluster cores over time, are likely to 
be rather complex. Previous studies showed that these 
sbSP and cSFR models do not fully explain the evolu- 
tion of the colors of cluster galaxies. However, they rep- 
resent valuable parametrizations of the observed CMRs, 
allowing for comparisons to previous work and direct, if 



approximate, estimation of ETC stellar population ages 
without the need for a complete specification of their var- 
ied histories. Therefore the above models will be applied 
as the initial framework for analyzing our results. 

5. RESULTS 
5.1. Color and Scatter Evolution 

We compile the measured colors and scatters of cluster 
CMRs in Table [3l and plot them in Figure [5] We also 
plot the associated sbSP and cSFR models, where for the 
latter we have utilized the parameter Zond rather than 
a weighted formation redshift (zj) because for a given 
model {zf) depends on the observation redshift. Thus 
the cSFR models plotted correspond to fixed histories, 
where star formation continues uniformly in time and 
ceases at Zqj^^ . 

In the top two rows of Figure \5\ we see that the rest- 
frame U — V colors of cluster RSGs are nearly constant 
across the redshift range 1 to 1.5, in contrast to any of 
the model stellar populations that redden monotonically 
with time. The color scatters in the bottom row of Fig- 
ure [5] reflect the same physical scenario: the rest-frame 
cr(C/ — V) is nearly constant (or slightly declining) with 
increasing redshift, while the scatter increases monoton- 
ically in cSFR models. 

We checked several other stellar p o pulation models 
jMa rastonI [20051: iConrov et all [20091: iConrov fc Guiii] 
12010) and parameter choices (metallicity, IMF, SFH) 
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Table 3 

Cluster red sequence color zeropoints and color scatters. 



Observed bands Rest-frame U-V 

Name z Vis. ETGs Quant. ETGs All RSGs Vis. ETGs Quant. ETGs All RSGs 



Color Zeropoint, Co 



ISCS 


J1429. 2+3357 


1.059 


1.92 




0.03 


1.93 




0.03 


1.93 ± 


0.03 


1.03 




0.03 


1.04 




0.03 


1.04 




0.03 


ISCS 


J1432.4- 


h3332 


1.112 


2.33 


± 


0.05 


2.32 


± 


0.06 


2.33 ± 


0.06 


1.05 




0.04 


1.04 




0.05 


1.05 




0.04 


ISCS 


J1426.1- 


h3403 


1.136 


2.04 


± 


0.05 


2.00 


± 


0.04 


2.00 ± 


0.04 


1.02 


± 


0.05 


0.99 


± 


0.04 


0.99 


± 


0.03 


ISCS 


J1426.5- 


h3339 


1.163 


1.95 


± 


0.07 


1.85 


± 


0.10 


1.87 ± 


0.06 


0.91 


± 


0.06 


0.82 


± 


0.08 


0.84 


± 


0.05 


ISCS 


J1434.5- 


h3427 


1.243 


2.35 


± 


0.05 


2.27 


± 


0.04 


2.34 ± 


0.04 


0.94 


± 


0.04 


0.88 


± 


0.03 


0.93 


± 


0.03 


ISCS 


J1429.3+3437 


1.262 


1.56 


± 


0.02 


1.53 


± 


0.02 


1.54 ± 


0.02 


1.05 


± 


0.02 


1.01 


± 


0.02 


1.02 


± 


0.02 


ISCS 


J1432.6- 


h3436 


1.349 


1.79 


± 


0.03 


1.81 


± 


0.02 


1.79 ± 


0.02 


1.12 


± 


0.02 


1.14 


± 


0.02 


1.12 


± 


0.02 


ISCS 


J1433.8- 


h3325 


1.369 


2.05 


± 


0.04 


2.01 


± 


0.05 


1.91 ± 


0.06 


1.33 


± 


0.04 


1.29 


± 


0.05 


1.20 


± 


0.06 


ISCS 


J1434.7- 


h3519 


1.372 


1.68 


± 


0.10 


1.71 


± 


0.08 


1.62 ± 


0.05 


0.97 


± 


0.10 


1.00 


± 


0.07 


0.91 


± 


0.05 


ISCS 


J1438.1- 


h3414 


1.413 


1.89 


± 


0.04 


1.85 


± 


0.03 


1.78 ± 


0.03 


1.11 


± 


0.03 


1.07 


± 


0.03 


1.00 


± 


0.03 


ISCS 


J1432.4- 


h3250 


1.487 


2.53 


± 


0.05 


2.49 


± 


0.06 


2.46 ± 


0.05 


1.01 


± 


0.04 


0.97 


± 


0.04 


0.95 


± 


0.03 



Color Scatter, (Tint 



ISCS 


J1429. 2+3357 


1.059 


0.08 




0.02 


0.12 




0.03 


0.13 ± 


0.02 


0.08 


± 


0.02 


0.11 


± 


0.03 


0.12 


± 


0.02 


ISCS 


J1432.4+3332 


1.112 


0.12 


± 


0.02 


0.20 


± 


0.07 


0.19 ± 


0.05 


0.09 


± 


0.02 


0.15 


± 


0.05 


0.14 


± 


0.03 


ISCS 


J1426.1+3403 


1.136 


0.18 


± 


0.03 


0.17 


± 


0.02 


0.18 ± 


0.02 


0.15 


± 


0.02 


0.14 


± 


0.02 


0.16 


± 


0.01 


ISCS 


J1426.5+3339 


1.163 


0.25 


± 


0.03 


0.30 


± 


0.06 


0.35 ± 


0.03 


0.21 


± 


0.03 


0.26 


± 


0.05 


0.29 


± 


0.02 


ISCS 


J1434.5+3427 


1.243 


0.17 


± 


0.02 


0.17 


± 


0.02 


0.19 ± 


0.02 


0.13 


± 


0.02 


0.13 


± 


0.01 


0.14 


± 


0.02 


ISCS 


J1429.3+3437 


1.262 


0.09 


± 


0.01 


0.09 


± 


0.01 


0.10 ± 


0.01 


0.09 


± 


0.01 


0.09 


± 


0.01 


0.10 


± 


0.01 


ISCS 


J1432.6+3436 


1.349 


0.05 


± 


0.01 


0.06 


± 


0.01 


0.07 ± 


0.01 


0.05 


± 


0.01 


0.06 


± 


0.01 


0.07 


± 


0.01 


ISCS 


J1433.8+3325 


1.369 


0.08 


± 


0.02 


0.06 


± 


0.04 


0.24 ± 


0.03 


0.08 


± 


0.02 


0.06 


± 


0.03 


0.22 


± 


0.03 


ISCS 


J1434.7+3519 


1.372 


0.28 


± 


0.06 


0.29 


± 


0.04 


0.23 ± 


0.03 


0.27 


± 


0.05 


0.28 


± 


0.04 


0.22 


± 


0.03 


ISCS 


J1438.1+3414 


1.413 


0.13 


± 


0.02 


0.12 


± 


0.02 


0.17 ± 


0.02 


0.12 


± 


0.02 


0.11 


± 


0.01 


0.16 


± 


0.02 


ISCS 


J1432.4+3250 


1.487 


0.15 


± 


0.03 


0.13 


± 


0.04 


0.21 ± 


0.02 


0.10 


± 


0.02 


0.10 


± 


0.03 


0.15 


± 


0.02 



and found similar results. Models with ~ OAZq can 
reproduce the flat {U — V)o ~ 1 trend at these redshifts, 
but cannot account for a flat or declining color scatter 
trend, nor for the colors of z < 1 giant cluster galaxies 
(|Aragon-Sa lamanca" et al.lll993l : [Stanford ct al. 19981). 

Several clusters are outliers to the median trends ob- 
served in these figures, and will be discussed in more 
detail in Sections 15.41 and El These general trends imply 
that we are not approaching the primary epoch of star 
formation in the centers of these clusters, and/or that 
their formation histories were not as simple as those as- 
sumed by the simplest models. Colors of the CMRs are 
redder than, and the color scatters smaller than the sim- 
plest sbSP or cSFR estimates in which the CMRs get 
bluer and wider as one approaches their primary forma- 
tion epoch. 



The change in U-V color owing to this effect is ap- 
proximately 0.1 magnitudes, and is therefore insufficient 
to explain the entirety of the deviation we see from any 
one sbSP in Figure [5l leaving the colors and scatters 
roughly constant across the redshift range. In this se- 
ries of figures (Figures [SI [51 [3 and [S]) we plot the en- 
tire RSG sample as unfilled black squares to demon- 
strate the (typically inconsequential) difference from the 
morphologically-selected ETC sample. We conclude that 
the cluster ETG samples are growing across the entire 
range 1 < z < 1.5 by the incorporation of younger RSGs 
that are missed in the samples at higher redshift, not only 
because they were star forming (and hence still growing) 
or disk-like in the centers at z > 1.3, but also because 
RSGs/ETGs that grow or enter the central region have 
younger stars than those previously present. 



5.2. Effect of Morphology Selection 

We demonstrate the effect of the morphology cut on 
both colors and scatters in Figure [HI At z > 1.3, the col- 
ors of morphologically-selected ETGs are slightly redder 
than the full red sequences at the same redshift. This 
might refl ect the canonical in-situ proge nitor bias as de- 
scribed in Ivan Dokkum fc Franxl (|2001| ) , where galaxies 
drop out of the ETG sample at higher redshift because 
they are continually being transformed and as such are 
not necessarily identified as ETGs in a cluster's progen- 
itor. Progenitor bias, or a similar effect, may provide 
a natural explanation for the seemingly unintuitive way 
that the red sequences remain narrow and red, a scenario 
we consider further in Section [6l The rate of transforma- 
tions inside the cluster center appears to have slowed by 
z < 1.3, lending additional evidence that above z w 1.3 
we are beginning to see the late stages of evolution by 
the giant cluster galaxies. 



5.3. Inferred (Simple) Star Formation Epochs 

Although the simple models of Section 31 may not re- 
flect the evolving stellar populations of the ISCS clusters, 
we can nevertheless gain some insight by fitting the in- 
dividual formation redshift for each cluster. Using the 
sbSP and cSFR models, we interpret the derived ages as 
the average time elapsed since the last major epoch of 
star formation in the galaxies currently making up the 
red sequence. 

The trends of H5.1\ imply a tight positive correlation 
between a cluster's observed redshift and its inferred star 
formation epoch. Below z « 1.3, these cluster RSGs 
resemble a single burst with {zj) « 2.0 (or Zend ~ 1-5), 
while higher-redshift clusters have {zf) > 3.0 (zcnd }^ 
2.0). We plot these fitted values in Figure [3 Error bars 
are calculated by fitting the endpoints of the one-sigma 
color and scatter ranges in the same manner. These data 
are compiled in Tabled In the lower-redshift suhsample, 
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Figure 5. Rest frame colors and scatters of our clusters' CMRs, 
with associated models. The top, middle, and bottom figures com- 
pare our measurements against the sbSP color, cSFR color, and 
cSFR color scatter models, respectively (Section |4}. Here we em- 
ploy visual morphological classification for the solid black diamonds 
and show the results for the entire RSG subsample as the open 
squares. For clarity, in this and later figures we only plot error 
bars for the morphologically-selected sample (black diamonds) and 
not the entire RSG subs amp le — they are similar in magnitude and 
both are listed in Tables IbIsI We plot a linear fit to the data in each 
panel as dotted black lines. There is a clear disagreement between 
the continued evolution of any single sbSP or cSF'R model and the 
observed CMR trends. 

no clusters exist that are the evolutionary continuation 
of clusters at higher z along the cSFR or sbSP tracks. 

Figure [8] and Table [5] show the same information, but 
converted into stellar ages, which we interpret as the av- 
erage time since the last major star-formation epoch for 
galaxies on the red sequence. These trends are qualita- 
tively similar regardless of the property-model combina- 
tion applied, sbSP color, cSFR color, or cSFR scatter. 
The scatter measurements in the bottom panel of Fig- 
ures [7] and [8] show the same general trend with redshift 
as the median color measurements of the top two pan- 
els: clusters in this sample at higher redshift indicate 
an earlier formation epoch. Therefore the time elapsed 
since the last major star forming eposide appears roughly 
constant, yielding a luminosity- weighted average age of 
« 2.3 ± 0.6 Gyr, consistent with the times cales inferred 
from rest-frame optical c olor scatters by IMevers et al.l 
(|2012D using the models of Ivan Dokkum fc Franxl (|200lD 
for seven of these clusters. 

Numerically, the best-fit ages show a large scatter, 
~ 0.5 to 1 Gyr, between when they are based on the scat- 
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Figure 6. Demonstration of the effect of our morphological se- 
lection, in terms of colors (top) and scatters (bottom), llere we 
plot the change in measured CMR quantities when we turn off our 
morphology cut. In the top panel, a negative number occurs when 
the CMR of ETGs is on average redder than CMR of all RSGs in 
the cluster. In the bottom panel, a positive number occurs when 
the CMR of ETGs is on average narrower than the CMR of all 
RSGs in the cluster. Although weakly, the data imply the exis- 
tence of bluer /younger non-ETGs on the red sequence in clusters 
above z ~ 1.3. One interpretation of this trend is as evidence of on- 
going morphological transformations of relatively young late-type 
galaxies into cluster ellipticals. However, these in-situ color dif- 
ferences do not account for all evolution in the median formation 
epoch seen in Figure [5] Therefore, a significant fraction of ETC 
luminosity at 2 ^ 1 was assembled from outside the cluster cores 
(or at least, outside the ETGs' progenitors) since z ~ 1.5, and 
composed of relatively younger stellar populations. 

ter and when based on median color. It may be that the 
simple star formation histories we assumed were inade- 
quate. For example, at a fixed red sequence color (and 
hence "age"), the color scatter between red sequence 
members can be made arbitrarily small by increasing the 
extent to which their star formation histories are corre- 
lated. Thus one possibility is that the formation histories 
of galaxies in some clusters are more uniform than others. 

5.4. Qualitative Discussion of CMR Assumption 

Several additional features pertaining to cluster galaxy 
evolution are seen in these data. Specifically, since the 
cluster selection technique does not depend directly on 
the presence of a canonical red sequence, these CMRs 
have a wide variety of properties. 

In nearly all cases, the cluster red sequence spans a 
factor of at least ^ 10 in brightness, indicating that 
cluster ETGs spanning the range L Ki 0.3 to 3 L* are 
present by z 1.5. On the other hand, we mea- 
sure several clusters with very tenuous or wide RSs. 
JSCS J1426.5+3339 (z = 1.163), JSCS J1434.7-I-3519 
(z = 1.372), JSCS J1433.8-F3325 (z = 1.369), and 
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Table 4 

Inferred CMR formation epochs. 



Name 



Vis. ETGs 



Quant. ETGs 



All RSGs 



Vis. ETGs 



Quant. ETGs 



All RSGs 



Prom sbSP Colors 



JSCS J1429. 2+3357 


1 


059 


2 


,„+0.14 
^'-'-0.14 


2 


12t« 


17 

18 


2 


15 


ISCS J1432.4+3332 


1 


112 


2 




2 


28t;] 


21 
27 


2 


38 


ISCS J1426. 1+3403 


1 


136 


2 


9/+0.28 
^+-0.24 


2 


05t;; 


20 
13 


2 


04 


ISCS J1426.5+3339 


1 


163 


1 


ocr+0.19 
°'J-().29 


1 


54t;; 


29 
04 


1 


55 


ISCS J1434.5+3427 


1 


243 


2 


if.+0.14 
l'J-0.12 


1 


89lE! 


13 
22 


2 


07 


ISCS J1429.3+3437 


1 


262 


2 


O1+0.15 
°1-0.19 


2 


50lg 


19 
15 


2 


62 


ISCS J1432.6+3436 


1 


349 


3 


c-7+0.21 
>"-(). 14 


3 


671° 


25 
13 


3 


55 


ISCS J1433.8+3325'' 


1 


369 




> 6 




> 6 




8 


86 


ISCS J1434.7+3519 


1 


372 


2 


c:/+0.S8 
°+-0.64 


2 


nl?, 


72 
39 


2 


27 


ISCS J1438.1+3414 


1 


413 


3 


7Q+0.23 
'^-0.19 


3 


561° 


15 
27 


2 


82 


ISCS J1432.4+3250 


1 


487 


3 


14+0-49 


2 


841° 


37 
22 


2 


73 



+0.14 
-0.15 
+ 0.15 
-0.28 
+0.18 
-0.12 
+0.23 
-0.02 
+0.11 
-0.11 
+0.17 
-0.17 
+0.12 
-0.10 
+0.91 
-5.03 
+0.21 
-0.39 
+0.30 
-0.19 
+0.20 
-0.16 



From cSFR Colors 



From cSFR Scatters 



ISCS J1429. 2+3357 


1 


059 


2 


l7^ 


-0.16 
-0.10 


2 


l9^ 


-0.21 
-0.15 


2 


2i: 


hO.19 
-0.11 


1 


oQ+0.09 


1 


+ 1-0.09 


1 


42: 


hO.lO 
-0.07 


ISCS J1432. 4+3332 


1 


112 


2 


44^ 


-0.25 
-0.24 


2 


33^ 


-0.34 
-0.20 


2 


46: 


F0.31 
-0.26 


1 


56+"-"' 
"^"-0.14 


1 


cn+0-20 

J"-0.12 


1 


57: 


1-0.19 
-0.15 


ISCS J1426. 1+3403 


1 


136 


2 


31^ 


-0.32 
-0.17 


2 


17^ 


-0.15 
-0.10 


2 


17: 


^0.13 
-0.09 


1 


J-'-o.ii 


1 


4-1+0.09 
+ 1-0.06 


1 


4i: 


h0.()8 
-0.05 


ISCS J1426.5+3339 


1 


163 


2 


01^ 


-0.17 
-0.10 


1 


88^ 


-0.11 
-0.08 


1 


89: 


^0.09 
-0.06 


1 


o., + ().09 
■'•^-0.06 


1 


9c:+0.07 
^■^-0.03 


1 


26: 


h0.()4 
-0.03 


ISCS J1434.5+3427 


1 


243 


2 


22^ 


-0.13 
-0.08 


2 


ll^ 


-0.07 
-0.08 


2 


2i: 


F0.08 
-0.08 


1 


-0+0.07 
48-0.05 


1 


40+0-03 

^^-0.04 


1 


46: 


ho.oe 

-0.05 


ISCS J1429. 3+3437 


1 


262 


2 


84^ 


-0.20 
-0.19 


2 


53^ 


-0.19 
-0.10 


2 


64: 


^0.18 
-0.13 


1 


R/1+0-13 
=4_o.ll 


1 


RS + 0-10 


1 


73: 


ho.io 

-0.07 


ISCS J1432. 6+3436 


1 


349 


3 


89^ 


-0.45 
-0.30 


4 


16^ 


-0.39 
-0.34 


3 


85: 


^0.30 
-0.22 


2 


^2 + 0.34 


2 


oq + 0.33 
°'^-0.26 


2 


59: 


hO.23 
-0.18 


ISCS J1433.8+3325=- 


1 


369 




> 


6 




> 


6 


8 


87: 


F0.23 
-4.50 




> 6 




> 6 


9 


22: 


hO.55 
-6.23 


ISCS J1434.7+3519 


1 


372 


2 


63^ 


-0.79 
-0.33 


2 


75^ 


-0.70 
-0.32 


2 


4o: 


F0.19 
-0.12 


1 


7f.+0.51 
'"-0.19 


1 


S/l + O-46 


1 


62: 


hO.ll 
-0.07 


ISCS J1438.1+3414 


1 


413 


3 


98j 


-0.54 
-0.39 


3 


52j 


-0.29 
-0.26 


2 


82: 


F0.30 
-0.13 


2 


6Q+0-42 
"^-0.30 


2 


O.+0.22 
•^+-0.17 


1 


91: 


h0.17 
-0.09 


ISCS J1432.4+3250 


1 


487 


3 


lfi+0.36 

^°-a.30 


2 


O7+0.34 

°'-a.i7 


2 


on+0.14 


2 


1 q + 0.26 

^■^-0.18 


1 


Qfi+0-21 
^"-0.12 


1 


qn+O-ll 



ISCS J1429.2+3357 


1 


059 


2 


01+0.17 

°l^-0.14 


2 


44+0-26 
++-0.23 


2 


qQ+0.18 

•^='-0.19 


1 


77+0.12 
' '-0.09 


1 ^54+0.17 
l-^+-0.12 


1 


CI +0.10 
■^1-0.10 


ISCS J1432.4+3332 


1 


112 


2 


7q+0.12 
'^-0.12 


2 


Q7+0.54 
^'-0.49 


2 


9-1+0.43 
^1-0.31 


1 


77+0.09 
' '-0.07 


1 47+0.32 
^•^'-0.28 


1 


40+0.25 

+•^-0.19 


ISCS J1426.1+3403 


1 


136 


2 


-,0+0.41 
!''-(). 17 


2 


97+0.33 
^'-0.21 


2 


-19+0.19 
1^^-0.10 


1 


QQ+0.24 
■^°-0.10 


1 47+0.19 

l-^'-0.14 


1 


07+0.12 
'"-0.05 


ISCS J1426.5+3339^ 


1 


163 


1 


70+0.22 
'■^-0.55 




< 2 




< 2 


1 


9n+o.09 

^"-0.08 


< 1.2 




< 1.2 


ISCS J1434.5+3427 


1 


243 


2 


_^ + 0.33 
'"^-0.34 


2 


7t-+0.28 
"J-0.32 


2 


4^+0.26 
^1-0.12 


1 


7q+0-20 
'=^-0.20 


1 7q+0.19 
1- ' ^-0.19 


1 


5g + 0.16 

■^''-0.08 


ISCS J1429.3+3437 


1 


262 


3 


40+0-08 

^"-0.13 


3 


QQ + 0.10 

•^^-0.15 


3 


24+0.14 
^^-0.26 


2 


90+0.07 
^■^-0.08 


9 99+0.08 
^•^^-0.10 


2 


-19 + 0.10 
1^-0.18 


ISCS J1432.6+3436 


1 


349 


4 


on+0.12 

•J^-cii 


4 


00+0.09 
•^•^-0.15 


3 


^'-0.13 


2 


70+0.11 
'"-0.09 


9 7r,+0.09 
^•'"-0.09 


2 


fi7+0.09 
"'-0.09 


ISCS J1433. 8+3325 


1 


369 


3 


qq+0.16 
^^-0.43 


4 


10+0.37 
l°-0.34 


2 


oq+0.17 
^^-0.25 


2 


00+0.12 
O°-0.31 


9 01+0.33 
''•°^-0.26 


1 


40+0.07 
+'^-0.35 


ISCS J1434.7+3519'' 


1 


372 




< 2 




< 2 


2 


-19+0.19 
-■-^-0.27 




< 1.5 


< 1.5 


1 


47+0.09 
+ '-0.36 


ISCS J1438. 1+3414 


1 


413 


3 


on+0.38 
^^^-0.39 


3 


07+0.42 
•^'-0.27 


2 


fii+0.11 
"1-0.12 


2 


-1 0+0.25 

l'^-0.23 
10+0. 58 
+•^-0.39 


9 9C.+0.29 
^■^■^-0.20 


1 


70+0.07 
'"-0.07 


ISCS J1432. 4+3250 


1 


487 


3 


fiQ+0.77 


3 


0-1 +0.98 
°l-().73 


2 


70+0.14 
'="-0.14 


2 


9 t;o+0.77 
Z.OO_o 4g 


1 


oq+0.11 

°y-o.o8 



In some cases, fitting resulted in seemingly unphysical parameters — e.g. galaxies formed earlier than our code 
permitted (~ 5 X 10^ years) or have a negative age — in these cases we give a one-sided estimate. 



ISCS J1432.4+3250 {z = 1.487) all contain an over- 
density of red galaxies that are near the expected red 
sequence color but do not define the canonical locus in 
color-magnitude space consistent with a single CMR of 
old galaxies. 

Moreover, clusters ISCS J1429.2+3357 {z = 1.059), 
ISCS J1432.6-I-3436 {z = 1.349), and ISCS J1434.5-I-3427 
(z ~ 1.243) appear to have red sequences where the 
red sequence is deficient in ETGs at L < 0.3 L* . 
Furthermore, ISCS J1434.5+3427 (z = 1.243) and 
ISCS J1438.1+3414 (z = 1.413) have 5 to 10 RSGs at 
L < L* that are bluer by up to ~ 0.4 magnitudes than 
the most apparent CMR, comprising f« 30% of the red 
sequence galaxies. We also notice a large number of 
early-type galaxies that are as bright and brighter than 



the red sequence we identified for ISCS J1438. 1-1-3414 
{z = 1.413), and that are ~ 1 magnitude bluer. How- 
ever, they are likely interlopers based on their projected 
spatial offset (> 0.5 Mpc) from the cluster center and 
association with a known foreground structure. 

Taken together with the quantitative results described 
above, these qualitative properties outline a picture 
whereby some massive cluster ETGs (ellipticals and SOs) 
are in place in many clusters, but perhaps not all, by 
z Si 1.5, with stellar populations formed at (zf) > 3. At 
z ~ 1, the stellar populations of the ETGs in this sample 
appear more recently formed, on average by w 1 Gyr, 
having {zf) ^ 2. 

6. DISCUSSION 
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Table 5 

Inferred CMR ages, in Gyr. 



Name 



Vis. ETGs Quant. ETGs 



All RSGs 









Stellar Age, 


^obs 


- (*/)' 


from sbSP Colors 




1 


059 


2 


30lg 


21 

23 


2 


qn+D.23 

•^■^-0.29 


2 


qO + 0.19 

■^"-0.23 


ISCS J1432.4+3332 


1 


112 


2 


66l« 


17 
34 


2 


■^'-0.39 


2 


f5Q+0.17 
"3 ="-0.37 


ISCS J1426.1+3403 


1 


136 


2 


44lg 


34 
36 


2 


J^O-0.21 


2 


-Tc+0.26 
^"3-0.19 


ISCS J1426.5+3339 


1 


163 


1 


731° 


33 
61 


1 


'-'='-0.11 


1 


in+0.49 
^'J-0.06 


ISCS J1434.5+3427 


1 


243 


1 


891° 


19 
19 


1 


°"^-0.43 


1 


OO+0.16 
°'J-0.17 


ISCS J1429.3+3437 


1 


262 


2 


661° 


13 
18 


2 


j-Ci 9n 
■^^-0.18 


2 


47+0.17 
-0.19 


ISCS J1432.6+3436 


1 


349 


2 


95lS 


11 
09 


3 


m +0.13 
'-'^-0.08 


2 


04+0.07 
y^-0.06 
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In some cases, fitting resulted in seemingly unpliysical parameters — e.g. galaxies 
formed earlier than our code permitted ('^ 5 X 10* years) or have a negative age — 
in these cases we give a one-sided estimate. 



We analyzed the rest-frame optical CMRs of red 
sequence galaxies in 11 spectroscopically confirmed, 
infrared-selected massive galaxy clusters at 1.0 < z < 
1.5. Clustering (B rodwi n ct al. 200J3), exten ded X-ray 
emiss ion ((Brodwin et al.H 2 01lf). weak lensing (iJee et al.l 
120111 ). and stellar masses ifEisenhardt et al.|[2008[) sup- 
port the view that these are massive (M > lO^^M©) 
clusters and are the precursors to present-day massive 
(M ^ IO^^Mq) clusters. These observations thus pro- 
vide a valuable window into the evolution of galaxies 
within the most massive halos in the universe by re- 
vealing the time during which most of their stellar mass 
formed and assembled from smaller pieces. 

6.1. Star Formation Histories of Cluster Ellipticals 

The present data (Figure [7]) indicate that massive, red 
ETGs existed in the centers of many rich clusters by at 
least z = 1.5. Other observations suggest this is true 
at even higher redshifts, both in massive, mature clus- 



ters ((Stanford et al.ll2012D and assembl ing structures or 
"proto-clusters" ('e.g.. lKurk et al.ir2009D . The red colors 
of these galaxies at z « 1.5 indicate that their luminosity- 
weighted age is « 2 to 3 Gyr (a formation epoch of 
z 4), and/or that the formation of stars then present 
ceased by ~ 1 Gyr prior to their observed state (z ^ 2.5). 

This story is rem iniscent of a "monolithic collapse" 
()Eggen et al.l 119621 ) . whereby the stars in these galax- 
ies formed in a single short burst at a much earlier 
time. Massive cluster ETGs at z < 1 evolve in a 
manner consistent with thi s passiv ely-evolving model 
(iBower et al.l[T99llEllis et al . 1997; S tanford et al.lll998l : 
iMancone et al.l 120101) . although their rest-frame optical 
colors are consistent with a wide variety of potential for- 
mation epochs (e.g., the first panel of Figure [5]) . 

The analysis presented here connects the low-redshift 
cluster galaxy population to their pasts by quantifying 
the stellar populations in cluster galaxies over a span of 
1.5 Gyr following the end of the peak star formation ac- 
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Figure 7. Inferred formation redshifts, corresponding to the same 
rows as Figure [5] The values in the top, middle, and bottom fig- 
ures are inferred from our sbSP color, cSFR color, and cSFR color 
scatter models, respectively (Section |4]l. Here we employ visual 
morphological classification for the solid black diamonds and show 
the results for the entire RSG subsample as the open squares. Clus- 
ter ISCS J1433. 8+3325 (z = 1.369) does not appear in the top two 
panels (nor in the visual ETG sample in the bottom panel) be- 
cause its CMR color indicates Zf > 6 (zond > 4). In all cases, 
the inferred red sequence formation epoch is later in lower-redshift 
clusters, suggesting a combination of ongoing accretion of younger 
stars by the massive cluster galaxies during this epoch and/or for- 
mation of younger clus ter ETGs. This extends the results found 
by, e.g. JMei et al.l l|2009 ) and Jaffc ct. al. {2011i) that the formation 
redshift inferred by SSPs increases with redshift. 

tivity in the universe at z ^ 2. The measurements of 
Section [5] suggest that the characteristic ages of the total 
steUar populations in massive RSGs are roughly the same 
across this cluster sample at 1 < 2; < 1.5 (FigureH]). This 
manifests as an inferred star formation epoch Zf that con- 
tinues to increase with the observed cluster r edshift, as 
found prev iously in samples p resented by, e.g.. lMei et alj 
(I2OO9I) and lJaffe et all mm . 

This is explained naturally if they experience signif- 
icant stellar mass assembly during this epoch, a sce- 
nario we discuss in Sections 16 .4H6 . 71 below. Alternatively, 
these observations could reflect drastic changes during 
this epoch in the nature of the galaxies/overdensities se- 
lected for this study, a possibility we consider in Sec- 
tions E^lJl 



Figure 8. Inferred ages, corresponding to the same rows as Fig- 
ure [5] The values in the top, middle, and bottom figures are in- 
ferred from our sbSP color, cSFR color, and cSFR color scatter 
models, respectively (Section |4]l. Here we employ visual morpho- 
logical classification for the solid black diamonds and show the 
results for the entire RSG subsample as the open squares. Clus- 
ter ISCS J1433. 8+3325 (z = 1.369) does not appear in the top two 
panels (nor in the visual ETG sample in the bottom panel) because 
its CMR color indicates an age > 4 Gyr. 



6.2. The Impact of Cluster Mass 

One factor that might confuse our analysis is if the 
nature of the clusters in the ISCS sample changes sig- 
nificantly over the redshift range we consider, or if they 
are fundamentally different structures than the progen- 
itors of today's massive clusters, as our interpretation 
assumes . The correlation fun ction of ISCS cluster can- 
didates (jBrodwin et al.ll2007| ) implies that they are pre- 
dominately M > structures and likely the pre- 
cursors to today's massive (~ lO^^Af©) clusters. Fur- 
thermore, for the spectroscopically confirmed subsample 
considered h erein, clusters with s ufficiently-deep X-ray 
observations (jBrodwin et al.ll20TT[ ) are detected at a sig- 
nificance consis tent with their being > IO^'^Mq objects 



Stanford et"aIll2012D . Finally, shear measure- 
Jee et al.l for all of these clusters with 



(see also, 
ments by 

very deep imaging from the HST Cl uster Supernova Sur- 
vey (PI Perlmutter, GO-10496; see lDawson et al.l 120091 ) 
yield individual (i.e., not stacked) weak-lensing masses 
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comfortably above IO^'^Mq. 

This does not preclude our spectroscopic subsample 
from containing a small number of structures that are 
somehow different than the remainder of the sample 
(for instance see JSCS J1433.8+3325 (z = 1.369) and 
JSCS J1429. 2+3425 {z = 1.162)). Furthermore, the ex- 
pected rapid assembly of cluster-mass halos at these red- 
shifts implies significant changes to the cluster galaxy 
population resulting from the incorporation of external 
groups via mergers, changes that might occur stochas- 
tically and lead to the outliers we observe in particular 
cases. Therefore the stellar-mass selection employed for 
this sample is effective at selecting cluster-size objects 
whose galaxy populations are not required to be fully 
evolved, permitting the analysis of a more diverse set 
of evolutionary paths in massive clusters. The fact that 
these clusters, which have similar stellar masses owing 
to their selection, exhibit similar colors implies that they 
had similar rates of star formation in their recent pasts. 
This matches our conclusion that massive cluster galaxies 
star formation over an extended period at 1 < z < 1.5. 

6.3. The Effect of Reddening 

Significant dust attenuation in cluster galaxies would 
bias our measurements toward a redder color, although 
we have no reason to believe that these ETGs contain an 
unusual amount of dust. If we posit that galaxies follow 
an evolutionary path from star-forming disks, through 
quenching, and ultimately to red-and-dead galaxies, then 
it is not clear that ETGs at higher redshift will have sig- 
nificantly mo r e dus t obscuration than those at z ^ 1. 
iMevers et al.l ()2012|) analyzed the ETGs in a set of clus- 
ters that includes seven members of the ISCS, and found 
that these galaxies, some of which are on the red se- 
quences in the present work, likely experience little dust 
reddening: E{B -V)< 0.06. 

If cluster ETGs at z > 1.3 are more highly ob- 
scured than their lower-redshift counterparts, but re- 
flect the same general star formation histories (suggest- 
ing zf 2), then cluster galaxies wo uld appear red- 
der (Eisenhardt et al.l l2008l and fainter (IMancone et al.l 
12010 ) . Line emission may also bias the present measure- 
ments to the red, since the WFC3/F160W filter includes 
the Ha line for clusters at 1.14 < z < 1.55. However, 
WFC3/G141 spectra (Zeimann et al. in prep.) indicate 
weak Ha emission in most of these RSGs. Furthermore, 
if either of these effects cause the red colors, we would 
expect significant variation between objects, and so the 
red sequence color scatters would systematically increase 
with redshift, contrary t o the constant or decreasing scat- 
ters we found in Section [5731 

This does not mean that dust plays no role at these red- 
shifts. At least one cluster's CMD (ISCS J1434.5-I-3427 
(z = 1. 243)) contains a spectroscopically confirmed 
member (jWagg et al.ll2012D satisfying the dust-obscured 
galaxy (DOG) criteria of iDev et afl ()2008[ ). In addi- 
tion, some objects (spectroscopic members or otherwise) 
are very red members of the clusters' red sequences. If 
this reflects contamination of the red sequences by dusty 
galaxies then our colors and/or color scatters may be bi- 
ased high in some cases, an effect we expect to get larger 
at higher redshifts owing to the increasing presence of 
dusty starbursts. This may also contribute to several of 
the outliers mentioned in Section [5.41 



Importantly, if dust or line emission do play a role 
in contaminating the red sequences or reddening clus- 
ter ETGs at z ^ 1.5, then it can only strengthen our 
conclusions that we are beginning to witness a period 
during which significant star formation activity recently 
occurred in massive central cluster galaxies. 

6.4. Luminosity Evolution 

Star formation over an extended period implies that 
cluster red sequences should be growing in stellar mass, 
and this implies evolution in the galaxy luminosity func- 
tion shape or normalization, scenarios that have been ex- 
plored in the ISCS by exploitin g the rest-frame n ear-IR 
data obtained with Spitzer. Mancone et al.l (|2010[ ) found 
that above z ^ 1.3, cluster galaxies are « 0.5 magnitudes 
fainter at 3.6 and 4.5/Ltm than one would expect from the 
extrapolation of their lower-redshift evolution, positing 
significant ongoing assembly at this epoch, leading to a 
rapid increase in luminosity with time. This conclusion, 
as well as the one cons istent with the present work by 
[Eisenhardt et al.l ([20081 ). is based on a sample of clus- 
ter galaxies that does not apply a color or morphological 
selection, lending more evidence to the notion that stel- 
lar mass must be actively increasing in cluster member 
progenitors. Rather than the entirety of such evolution 
being caused by progenitor bias, we are likely witness- 
ing the effect of real growth in the stellar mass of some 
cluster galaxies. 

The results b^ [Mancone et al.[ (|2010( ) and 

'Eisenhard t et al.l (|2008[ ) were based on photometric 
cluster members within 1.5 and 1.0 projected Mpc, 
respectively, of the cluster candidate's center. Like we 
posit for the present work, those studies may track 
a galaxy sample that misses stellar mass that formed 
outside of these regions, or in galaxies not bright enough 
to be selected. Thus we suggest the luminosity growth 
in massive cluster galaxies at these redshifts is due in 
part to star formation that occurred over an extended 
period, rather than arising entirely from the merging of 
equally old progenitors. 

6.5. Red Sequence Growth 

Then what is the star formation and assembly history 
of massive cluster galaxies? The data in Figure [5] are 
consistent with each cluster's ETGs forming at a differ- 
ent time, as the sbSP models with Zf — 2 to 6 roughly 
bracket the observed colors. We might reasonably expect 
that the star formation epoch of the massive galaxies in 
one cluster differs from another based on the assembly 
history or halo mass, and that differences in this history 
are small enough to have damped out by z < 1. 

However, cluster galaxies in the present sample at 
1.0 < z < 1.3 tend to be bluer than those implied by the 
continued passive evolution of the clusters at z > 1.3: 
any sbSP {zf > 3) that fits the high-redshift clusters 
is too red for the lower-redshift ones, suggesting that 
the colors of massive cluster galaxies might evolve in 
a manner altogether different from the commonly-used 
simple sbSP and time-uniform composite models consid- 
ered in Sec tion H) Such evolution has been o bserved for 
some time ([Beh et al.[[200l [Brown eTsIMOTh m surveys 
of massive field galaxies. Additional significant growth 
of the massive red galaxy population occurs at z > 1 
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(iKriek et al.ll2008Hvan Dokkum et aIll2008HTavlor et alj 
l2009HWhitaker et al.ll2010D . and this may be likely true 
in all environments. We therefore consider the possibil- 
ity that our cluster sample selects, at least in part, a 
continuous evolutionary sequence that experiences star 
formation over an extended period leading to an increase 
in stellar mass on the red sequence. 

Importantly, for the ETGs to reside on the red 
sequence and maintain a nearly-constant luminosity- 
weighted age as we observe, the red sequence must ac- 
crete newer stars in some fashion, making the CMR bluer 
than one which evolves passively. Moreover, the newer 
stars or galaxies must have ceased forming some time ago 
(> 100 Myr) in order to be considered RSGs, and must 
comprise a body of stellar mass that was experiencing 
an increasing star formation rate prior to its quenching. 
Otherwise the observed luminosity-weighted age would 
increase with time. Thus, these stars likely formed in 
gas-rich, rapidly star-forming objects that are not found 
in large numbers on the red sequence in the centers of 
their progenitor clusters. In this picture, the halting of 
star formation in massive cluster ETGs does not occur 
suddenly and absolutely in every galaxy of the final clus- 
ter, but continuously, during or after periods of signifi- 
cant gas accretion in objects that eventually make up the 
final, present-day ETGs. 

With this interpretation, the OMRs of these clus- 
ters may reflect in part the ongoing growth of galaxies 
in the cluster centers as they accrete or form younger 
stars not present at a previous time, as well as the 
transformation from star-forming to quiescent galaxies. 
The latter is the canonical in-situ progenitor bias of 
Ivan Dokkum fc Franxl ()2001t ). who described this effect 
at z < 1 for satellite galaxies, which exhibit remarkable 
constancy in their optical colors and color scatters over 
time. We now see this same result in massive cluster 
ETGs (Figure [5]), and we showed in Section that the 
transformations of late-type RSGs into ETGs are not 
sufficient to account for the entirety of the evolution we 
observe in the CMR formation epochs. Therefore an- 
other form of progenitor bias likely contributes, in which 
the red galaxy samples in the cluster centers gain stellar 
mass over time from outside the RSG sample in order to 
achieve roughly constant colors and scatters. 

Contrast this result with a mass-selected field galaxy 
samp le at compara ble rcdshifts such as that considered 
by W hitaker et al.l (2010), who found that the color scat- 
ter of field RSGs increases rapidly from z = ltoz=1.5. 
This likely owes to the way our sample of massive cluster 
galaxies were selected: we require them to exist in the 
center of a cluster, the bottom of a deep potential well 
that continuously collects significant amounts of stellar 
mass. We see only those galaxies and stars present there 
at the time we observe them, but may miss a significant 
amount of stellar mass that formed in massive galaxies 
outside the central ETGs at z < 2 to 3. 

6.6. Models of Ongoing Assembly 

We now seek an understanding of the impact of 
stellar mass assembly on the appearance of cluster 
cores over this redsh i ft ran ge. Previous work (e.g., 
Ivan Dokkum fc Franxl I2001D highlighted the effect of 
"progenitor bias", whereby galaxies identified as ETGs 
or RSGs at one redshift are not selected in compara- 



ble higher or lower redshift samples. Potential causes 
of this effect include the continued formation or growth 
of these galaxies, bursts of star formation in RSGs, 
and/or the transformation of star-forming disk galaxies 
into RSGs/ETGs via merging or other processes. There- 
fore we expect this and similar effects to be apparent in 
this sample. 

As galaxies form stars, possibly merge, and eventu- 
ally experience quenching as they fall into the center 
of the assembling halo, significant amounts of time may 
pass before they become incorporated into what we see 
as galaxies in the cluster center. Also, massive galaxies 
are thought to accrete a significant fraction of their stel- 
lar mass from their environment, thro ugh numerous mi- 
nor merger s or smooth accretion fe.g. .lKeres et al.| | 2005l : 
IWhite et a l. 2007; Dckcl ct al. 2009; Oser et al.ll2010t) . In 
this way massive galaxies experience assembly of struc- 
tures that alter their average stellar population with 
time. If the stellar ages of the new material are differ- 
ent enough from those already assembled, then arbitrary 
paths through the space of inferred ages (through a sim- 
ple proxy such as color) are possible, depending on the 
precise history of assembly. 

In any case, the salient result of this structure growth 
is that the cluster cores grow in stellar mass density, a 
generic process for which we seek an intuitive descrip- 
tion in terms of the quantities we measure here. A com- 
plete specification of the different star formation histo- 
ries in massive structures is beyond the scope this work, 
yet simple models can demonstrate the key character- 
istics expected of galaxies that evolve in this fashion. 
Fo r this purpose we app l y thre e such models inspired 
by Ivan Dokkum fc Franxl (|200lD that differ only in the 
time delay between when a galaxy's stars are formed and 
when they are incorporated into the ETC sample. Model 
galaxies that have formed within Ai of the observation 
time are masked from the corresponding CMR model, 
mimicking the time before a galaxy would be identified 
as an ETC or RSG in the center of a cluster. Thi s proce- 
dure has been used previously (for example. ' Jaf fe et al.l 
120111) to analyze the evolution of cluster galaxies using 
their red sequence colors and scatters. 

This delay, and its time evolution, could represent vari- 
ous physical scenarios, including galaxies that experience 
a small burst that removes them from the red sequence, 
galaxies that have not yet entered the center of the clus- 
ter, and galaxies that are not bright enough to be de- 
tected in the sample. We do not apply any color cuts 
on these model CMRs: all model galaxies are assumed 
to be on the red sequence. If a particular choice of At 
leads to very blue galaxies, then the CMR model would 
be relatively blue on average. We could test an infinite 
number of choices for At, but our data can only reason- 
ably constrain a small number of representative scenar- 
ios. Therefore we will test three such toy models in the 
hope of gaining some additional physical insight into the 
formation of massive cluster galaxies. 

For the first two models we choose to demonstrate a 
constant At of 1.0 and 2.3 Gyr, ro ughly consistent with 
the tdeiay ~ 2 Gyr value found by IMevers et al ' 
after utilizing models by Ivan Dokkum fc Fra: 
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to the optical color scatter in a subset of seven ISCS 
clusters. Our third model, "sudden infall" , corresponds 
to sudden mass growth by assuming that At is a func- 
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tion of redshift, with a minimum of 200 Myr occurring 
at z = 1.2 but steeply rising to 2.3 Gyr before and af- 
ter. Because galaxies that are younger than the existing 
RSGs are being added in large numbers to the CMR at 
1-2 ^ ^ 1-8, this corresponds to a model red sequence 
that becomes bluer with time, where the amount of color 
variation depends on the suddenness of the assembly, the 
minimum delay time, as well as the closeness of the in- 
fall time to the time of peak SFR in the galaxies that 
eventually make up the cluster. In this case, because we 
track an estimate of the universal SFR, this peak occurs 
at z ^ 2 to 3. 

To implement these models, we assume that each 
galaxy is formed in a single short burst and enters the 
cluster ETG sample in a first-made, first-in fashion ex- 
actly At after its formation. We simulate 10^ sbSP tracks 
that form at a rate such that the final red sequence expe- 
rienced the same star formation history in all three cases. 
To reduce the number of free parameters, we use mod- 
els with__ttiejjjuv^^ history estimated 
by iHopkins &: BeacomI ()2006f ) . For simplicity we simu- 
late only the galaxy color distribution and not the entire 
CMR. 

These models have the beneficial features that a) the 
colors of the galaxy populations in higher redshift clus- 
ters are redder than one would expect given a mono- 
lithic collapse scenario with a single formation redshift, 
b) these colors are consistent with an earlier inferred star- 
formation epoch without necessarily leaving a redder de- 
scendent, and c) the final model "red sequences" have 
identical stellar populations, regardless of the details of 
their assembly history. Therefore they only reflect the 
consequences of (highly idealized) structure growth and 
do not depend on specific assumptions for the physics of 
star formation quenching in massive halos. 

6.7. Colors and Scatters of Assembling Red Sequences 

We plot the color and color scatter evolution of these 
assembly- inspired models as the gray curves in Figure [9l 
In the top panel we summarize the intrinsic U-V colors 
of our sample and compare it to data from the litera- 
ture on clusters at 1.0 < z < 1.9. In this panel, colors 
are computed relative to the Coma cluster and denoted 
A{U - V)o. The color scatters of cluster ETGs (Fig- 
ure [9l bottom) agree largely with inferences made from 
their colors only, as seen in Figures (U [71 and [H 

As seen in previous work, the immediate result of the 
delay between the galaxies' formation and their settling 
on the red sequence is that cluster red sequences con- 
tinue to remain redder and narrower than any single sbSP 
or cSFR model implies, even while the red sequence is 
sparse or tenuous as observed in some cases fSection l5.4|) . 
presumably before they are fully formed. 

The feature that inferred ages appear more or less con- 
stant at z > 1 is a natural consequence of the stars enter- 
ing high-redshift ETG samples having been drawn from 
a history with an increasing star formation rate. This 
scenario wa s observed in di stant galaxi es (z > 3) by , 
for example. iLee et ahl (|2011[ ) and Papovi ch et al.l (|2011l ). 
suggesting that the galaxies in assembling clusters may 
be imprinted with this rising-SFR signature from their 
gas accretion epoch. In that case, for a given delay At, 
the dominant stellar population in the optical at any time 
resembles an sbSP formed a time At ago. Our third such 



model experiences a rapid decrease in the SF-assembly 
delay at z = 1.2, possibly mimicking the result of a dra- 
matic enhancement to star formation in galaxies falling 
into the cluster, such as a merger with a star-forming 
group or a temporary increase in smooth gas accretion. 
This demonstrates one plausible effect that can make as- 
sembling CMRs bluer (or wider) for a time (cf. cluster 
ISCS J1426.5-h3339 (z = 1.163)), which may be neces- 
sary or expected during this epoch. We do not claim that 
our data are yet able to constrain this type of model in 
detail, but these models do provide some evidence that 
such scenarios, and perhaps combinations thereof, might 
naturally occur in a sample like the ISCS. 

The color scatter measurements provide an additional 
constraint, because they reflect the extent to which the 
star formation episodes of the cluster ETGs (at fixed 
color) were correlated. We demonstrate specific variants 
of this in the bottom panel of Figure [SI the assembly- 
inspired models can yield the necessary slowly-varying 
(or increasing) color scatter as time progresses for clus- 
ters at z > 1. In these models, each "galaxy" is an 
uncorrelated sbSP formed at a random time such that 
the correct average SFH is borne out in a large sample. 
Therefore the correlation between model ETG star for- 
mation episodes reflects only the correlation implied by 
the underlying SFH (e.g., when it peaks, galaxies form- 
ing at that time will have correlated SFHs), and so any 
such scatter models could be considered as upper limits. 

The model tracks we show in this work are « 0.05 
magnitudes bluer than clusters at z < 1, reflecting either 
a possible systematic offset in the photometric systems 
between the model and sample (here we merely quote 
the numeric values from the cited works), that cluster 
galaxies experience an earlier peak in their star forma- 
tion than the SFH we assumed, or that massive cluster 
ETGs at low redshift have a slightly higher than solar 
metallicity. 

While order-unity departures from roughly solar metal- 
licity in cluster ETG stellar populations are inconsistent 
with their evolution at z < 1, the metallicity evolution in 
cluster ETG progenitors at z 1 is unclear. The color 
and scatter trends we have observed here, at a fixed for- 
mation epoch {zf ^ 3) using our sbSP or cSFR models, 
would imply that the luminosity- weighted average stellar 
metallicities decrease by a factor of ~ 5 during approxi- 
mately 1 Gyr at 1 < z < 1.5. This by itself is perhaps un- 
physical, but if it does contribute, it again implies stellar 
mass growth, in this idealized case owing to the incorpo- 
ration of relatively low-mass (low-metallicity) galaxies. 
Like our toy models at fixed metallicity but evolving av- 
erage age, this scenario would suggest significant ongo- 
ing star formation near cluster cores, and assembly into 
cluster cores, over an extended period at z ^ 1.5. In re- 
ality both effects might contribute to the observed trends 
at this epoch, when assembling clusters contain massive 
galaxies (high-Z, early z/), but grow their red sequences 
rapidly by incorp orating or growing low er-mass (lower-Z, 
late Zf) galaxies (jRaichoor et al.ir201lD . 

In this work, we have not separated the elliptical ver- 
sus SO galaxy samples, and therefore the apparent evo- 
lution of the ETG population may owe to an increasing 
population of SO galaxies (and not neces sarily ellipticals) 
with time, extending the results of, e.g.. iPoggianti et al.l 
(|2^00a) . Subsequent analyses will determine the evolution 
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Figure 9. Top: Evolution of the rest-frame U-V colors of z > 1 cluster galaxies, relative to Coma. Bottom: Rest-frame U-V scatters 
of cluster galaxies. We plot cluster data from the present work, and data compiled by Hilton et al. (2009)", including clusters from 
IBlakeslee eTaTI ||2003| ') and IMei et"all I l2006al lbl). We also include clusters CIG J021 8.3-05 10 (Papovich et al. 2010)^ and two new clusters 
from the IRAC Deep Cluster Survey'^, IDCS 1426+3508 {z = 1.75; [Stanford c t'^l20Tl ) and IDCS 1433+3306 {z = 1.89; Zcima nn et al.l 
120121 '). The curves are different toy models for the stellar populations in massive ETGs at the centers of clusters. The solid red, green, and 
blue curves are the cSFR models corresponding to Zend = 3, Zend = l-8i a^nd ^end = 1-2, which we described in Section |4] Wc show the 
simpl e alternatives of Section l6.6l in the gray curves: these models experience the same overall star formation histories (Hopkins & Bcaco^ 
120061) . but differ in the time delay between the formation of stars in a galaxy and the time at which the galaxy becomes apparent in the 
cluster red sequence. 
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of the cluster galaxies' morphologies. 

7. CONCLUSIONS 

Wc studied sensitive, high-resolution, rest-frame opti- 
cal imaging of the central regions in 13 spectroscopically 
confirmed infrared-selected galaxy clusters at 1 < ^; < 1.5 
that are plausible progenitors of today's massive clusters 
with M ~ 10^^ Mq. In 11 of these clusters we measured 
the CMR of red galaxies at a color consistent with the 
cluster's known redshift, and subsequently we analyzed 
the stellar populations implied by the location and dis- 
persion of the red sequence. We find the following: 

1. Red sequences of bright galaxies are present in the 
centers of many clusters by 2; = 1.5. However, the 
properties of their CMRs indicate large cluster-to- 
cluster variations in the formation histories of clus- 
ter galaxies observed at a given redshift. 

2. As a function of redshift from z = 1.5 to 1.0, the 
median rest- frame U — V color of cluster galaxies in 
this sample is nearly constant at ([/ — V)o ~ 1.0, 
and the scatter about the median color is nearly 
constant at CTint ~ 0.1 magnitudes. Some clusters 
deviate widely from these trends, and outliers may 
represent contamination by dust-obscured galaxies, 
the evolution of galaxies through the "green val- 
ley" , or red sequences that are not yet fully formed. 

3. The flatness of the color and scatter trends at 

1 < z < 1.5 requires star formation histories more 
diverse than a single short burst in the past. These 
average trends are plausibly in part the result of 
progenitor bias, where the progenitors of lower- 
redshift cluster ETGs are excluded by the sample 
selection at higher redshifts. 

4. Across this redshift range, these average colors and 
scatters both imply a roughly constant luminosity- 
weighted stellar age in the cluster RSGs, where 
the time since the bulk of stars formed is roughly 

2 to 2.5 Gyr, and the time since star formation 
halted is roughly 1 to 2 Gyr. 

5. A constant stellar age is a natural consequence of 
observing the center of an assembling halo whose 
final stellar population is experiencing ongoing or 
increasing star formation at the observed redshift. 
If the red sequence stars are nearly completely in 
place, they would reflect a small continuing star 
formation rate and their luminosity-weighted ages 
should increase with time. By contrast, at z ~ 
1.5 our observations imply that significant current 
star formation is occuring that will contribute to 
the stellar populations of central cluster galaxies 
by 2 ~ 1. 

6. Assuming the star formation of cluster galaxies fol- 
lows that of other systems, the color and scat- 
ter trends we find are consistent with a factor 
of two growth in the cluster red sequences from 
z = 1.5 to 1.0. 

7. These observations lend credence to the idea that 
the red sequence in the centers of clusters was grow- 
ing rapidly at z ~ 1.5, where this mass assem- 
bly was likely in the form of relatively younger, 



smaller, and/or metal-poorer galaxies that were ac- 
tively growing in their recent past. 
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